Spatio-temporal oscillations and rheochaos in a simple model of shear banding 
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We study a simple model of shear banding in which the flow-induced phase is destabilised by 
coupling between flow and microstructure (wormlike micellar length). By varying the strength of 
the instability and the applied shear rate, we find a rich variety of oscillatory and rheochaotic shear 
banded flows. At low shear and weak instability, the induced phase pulsates in width next to one 
wall of the flow cell. For stronger instability, single or multiple high shear pulses ricochet across 
the cell. At high shear rates we observe oscillating bands on either side of a defect. In some cases, 
multiple such defects exist and propagate across the cell to interact with each other. We discuss our 
results in the context of recent observations of oscillating and fluctuating shear banded flows. 
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Complex fluids commonly undergo flow instabilities 
and flow-induced phase transitions that result in spa- 
tially heterogeneous "shear banded" states. Classically 
studied systems include wormlike micellar surfactants 1]; 
dense lamellar onion |2| or micellar cubic 3] phases; and 
polymer solutions 0. Fluidity banding has also been 
reported in soft glassy materials Q such as colloidal sus- 
pensions 6] and simulated Lennard Jones particles Q. 
Experimentally, the basic observation is of two coexist- 
ing shear bands with differing viscosity and microstruc- 
ture (or fluidity). Theoretically, this is captured by in- 
voking multiple flow branches in the constitutive relation 
of shear stress vs. shear rate, £(7) j8l. Ho] . The system 
then separates into a steady state comprising two shear 
bands, each on its own flow branch (Fig.^Ji). 

However an accumulating body of data shows this ba- 
sic picture to be oversimplified: Many shear banding 
systems display oscillations or irregular fluctuations sug- 
gesting chaos in their bulk rheology, rheo-optics or ve- 
locimetry. Example systems include onion phases (at the 
layering transition or showing size oscillations fl2^: 
wormlike micelles (WMs) with a stress plateau in the 
flow curve [T^]; shear thinning WMs showing common 
stress banding (band normals in the flow-gradient direc- 
tion) Il4l IT^ : shear thickening WMs showin g; co mmon 
stress [l6j or common strain rate (vorticity) [T3 band- 
ing; and polymer solutions 0. In contrast, present mod- 
els predict steady banded states 0, 0, [HI • 

In the negligible Reynolds limit relevant to these vis- 
coelastic materials, the non linearity underlying this er- 
ratic response must arise in the rheological constitutive 
behaviour of the system |l8|. Temporal "rheochaos" has 
recently been studied in homogeneous models of director 
dynamics in sheared nematics ^?j; and of shear thick- 
ening systems with a single-branched flow curve [l8|- In 
many systems, however, spatial heterogeneity is likely to 
be a crucial ingredient of rheochaos. In this work, mo- 
tivated by the above experiments, we introduce the first 
model of spatio-temporal rheochaos in shear banding sys- 
tems with multi-branched constitutive curves. 



Homogeneous flow is unstable in any region of nega- 
tive constitutive slope, d^E/dj < 0. This is easily seen 
in models that take the "mechanical variables" (E, 7) as 
the relevant dynamic variables [2(J • The system can then 
separate into two shear bands, each on its own stable flow 
branch. (See Fig.^Ji, for shear thinning gradient-banding 
systems.) In more realistic models, the mechanical vari- 
ables are coupled to microstructural quantities such as di- 
rector orientation Q , polymeric concentration [2]| or mi- 
cellar length • This coupling can destabilize the rising 
high shear branch (causing, e.g., tumbling and wagging in 
nematics 23]). Here we construct a simple model with an 
unstable high shear branch (Fig. Q]3,c) and show that it 
has oscillatory and erratic (chaotic) shear banded states 
at imposed global shear rate. Although we have not yet 
proved the erratic states to be strictly chaotic, we refer 
to them as chaotic hereafter. Our model resembles other 
globally coupled reaction diffusion systems, with similar 
spatio-temporal behaviour |25|. Aradian and Cates are 
currently studying spatio-temporal flows of similar mod- 
els with a single unstable constitutive branch [2^ . 

We use only the minimal ingredients needed to capture 
the observed phenomena, so do not address the micro- 
scopies of any given system. However for concreteness 
we use the language of shear thinning WMs. There ex- 
ist numerous reports of apparently unattainable homoge- 
neous high shear rate branches in such systems, in which 
the sample flows erratically or is ejected from the cell 
(e.g. 1]). This is seldom discussed in detail, but some- 
times attributed to surface instability. Nonetheless, the 
possibility of bulk instability remains. Indeed, in some 
WMs the high shear band breaks into sub-bands [26|. 

For simplicity we consider just one microstructural 
variable, the mean micellar length n, and define our 
model by 
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FIG. 1: Intrinsic constitutive curves for differing degrees of 
coupling between flow and micellar length, a) Weak coupling, 
giving the standard coexistence of stable low and high shear 
bands (A and B); b) moderate coupling; c) strong coupling. 
Squares show Hopf bifurcations. The thin black lines delimit 
the periodic orbit of the local model at fixed E. 

Eqn. n gives the uniform shear stress E(£) in the non- 
inertial limit. It comprises a (generally no n- uniform) vis- 
coelastic micellar contribution cr(y,t) and a solvent con- 
tribution with viscosity n. The dynamics of a (Eqn. EJ) 
has a length-dependent relaxation time [2?| 

r(n)=r (-)\ (4) 
\ n oJ 

and a steady homogeneous state a — g[jr(n)] set by 



admitting a constitutive curve of negative slope. The 
spatial gradients DdyCr allow for a smooth interface of 

width / oc \[D between the shear bands. 

The micellar length n has its own relaxation time r n 
(Eqn. [3J), presumably related to the (unknown) under- 
lying rates of micellar scission and recombination. For 
example, scission could be enhanced by the agitation of 
shearing; or recombination aided by shear-induced align- 
ment of micellar ends [22[. We assume the former, taking 

= tt^- ^ 

Because this decrease in length only becomes important 
for typical shear rates 7 > l/r n , the strength of cou- 
pling between the mechanical quantities (E, cr, 7) and mi- 
crostructure n can be tuned by varying r n . 

Using this model we study flow between two parallel 
plates at y = 0, L with boundary condition c^E = 0, 
using units in which no = l,ro = 1 and L = 1. We set 
a = 1.2, (3 = 1.5, though our results are qualitatively 
robust to reasonable variations in these values. 

Homogeneous dynamics — The constitutive curves 
£(7) and domains of instability for homogeneous states 
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FIG. 2: Homogenous dynamics (D = 0) at fixed logE = 
—0.301 for different initial conditions. r n = 0.145. 



are shown in Fig. E] For small r n < 0.115 we find pure 
mechanical instability (Fig. QJi) that does not involve n. 
The constitutive curve is then (at any fixed E) an unsta- 
ble saddle |2£| (one unstable eigenvector). For larger r n , 
coupling to micellar length broadens this instability into 
the rising high shear branch, which is now an unstable 
focus [23 (two unstable eigenvectors): Fig. ^>,c. This 
instability terminates in a Hopf bifurcation |28| . 

We then studied the non-linear dynamics of the local 
(D = 0) model for fixed E, solving Eqns. (IlloP via a 
fourth order Runge-Kutte method [29(. This confirmed 
the stability properties of Fig. ^ states near an unsta- 
ble (stable) segment of the constitutive curve flow away 
from (towards) that segment (Fig.[2J). They also reveal a 
periodic orbit about the unstable high shear branch for 
stresses just below the Hopf bifurcation. We also used 
AUTO |3(| to trace the amplitude of the periodic orbits, 
Fig.^3,c. Periodic orbits are the most complex behaviour 
possible for the local model since it has only two degrees 
of freedom, d = 2. Chaos requires d > 3 |28j |. 

Spatially heterogeneous dynamics — We now turn to 
the non-local model, D / 0, focusing on the impli- 
cations of an unstable high shear branch (Fig. H^c). 
The dimensionality d is now effectively infinite, since 
each spatial point has its own value of n and 7. We 
solved the non-local equations using a Crank-Nicholson 
algorithm (checking our results with the Rosenbrock 
method) 29], with the constraint of fixed global strain 
rate, 7 = J* dyy(y,t). For small r n < 0.115 we find 
stable shear bands (Fig. QJi). In contrast for r n = 0.145 
(unstable high shear branch) we find spatio-temporal os- 
cillations and chaotic banding flows (Fig.|3J). 

Several regimes are evident. At low applied shear, 
7 = 1.5, a thin pulse of high shear ricochets back and 
forth across the cell (Fig. 03 top). A thin fluctuating high 
shear band, away from the rheometer wall, was seen ex- 
perimentally in Ref. At larger shear rates, we find 
two or more such pulses. For two pulses (not shown), 
we typically find a periodically repeating movie with the 
pulses alternately bouncing off each other (mid-cell) and 
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FIG. 3: Right: space time plots for shear rate evolution for 
r n = 0.145, D = 0.0016 with space y = - 1 left-right and 
time t — 120 — 140 bottom-top. Greyscale for the shear rate 
is shown in each case. Left: corresponding stress vs. time. 

the cell walls. Once three pulses are present, 7 = 7.0, 
periodicity gives way to chaotic behaviour (Fig. |3J). 

At still higher shear, 7 = 19.2, we find regular oscilla- 
tions of spatially extended bands pinned at a stationary 
defect. The local shear rates span both the low and high 
shear constitutive branches. Oscillating (vorticity) bands 
were seen experimentally in Ref. [j, [l7| . For the inter- 
mediate value 7 = 11.35 we find intermittency between 



FIG. 4: Space-time data of Fig-Hfor 7 = 1.5, 11.35, 23.0, 31.0 
(arrows under abscissa) and 130 < t < 140 in condensed form. 
Solid line - intrinsic constitutive curve; long dashed line - 
extrema of the periodic orbits of the local model; dotted line 
- maximum shear rate of the propagating high shear pulse 
described by Eqns. |H| 



patterns resembling those for 7 = 7.0 and 7 = 19.2. 

Finally for 7 = 23.0, 31.0 we find oscillating bands sep- 
arated by moving defects (Fig. [31 bottom right) with the 
flow now governed only by the high shear constitutive 
branch. In each band, the shear rate cycles round the 
periodic orbit of the local model (Fig. [3 and top data set 
of Fig. HI). 

For different values of r n we find a host of other in- 
teresting phenomena |3]|. For example, for weaker in- 
stability (r n = 0.13) at low applied shear rates we see a 
high shear band that can pulsate in width while adher- 
ing to the rheometer wall (Fig. (St), or meander about the 
cell (Fig.JHt)). The former behaviour resembles interfacial 
motion in seen WMs 0? ^3 an d onion phases . 

Finally we discuss in more detail the high shear pulse of 
Fig. |3| (top right). At times when the pulse is far from the 
wall, the stress E is constant (flat regions in Fig. top 
left). In this regime we transform to the pulse's comoving 
frame y = y — ct and eliminate a to get 
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in the 3d space (71,7,7'), with parameters E and c. 
(Prime denotes d/dy.) The low shear regions in Fig. [31 
(top right) lie on the low shear branch of the constitutive 
curve. Although this is a stable stationary flow branch 
(Fig. [IJ, in the space of (71,7,7') it is actually a fixed 
point (P) with one unstable eigenvector and two stable 
ones. As y increases through the high shear pulse, the 
solution's trajectory moves away from P along the un- 
stable direction and then back to P in the stable plane, 
completing a "homoclinic connection" 28]. The require- 
ment for this Id trajectory to correctly rejoin the 2d sta- 
ble manifold selects (for a given E) the observed prop- 
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FIG. 5: Shear rate density plots for r n = 0.13, D = 0.0016. 

agation speed, c. For other values of c, the trajectory 
abruptly crosses the 2d manifold, missing P. (This fol- 
lows by simple dimension counting Once selected, 
the solution's profile j(y) can (for a given value of D) be 
integrated over y to find the globally applied shear rate 7. 
By the reverse reasoning, we find a selected E and c for 
any 7 (and D). We used AUTO 30] to check that the 
homoclinic connections of Eqn. [5| indeed coincide with 
our numerical results. The locus of the maxima 7h(E) 
of these connections is marked in Fig. 21 The space-time 



points for 7 > 7h arise during wall collisions; correspond- 
ingly, their density is smaller than for 7 < 7h- 

To summarise, we have constructed a simple model 
of shear banding in which the high shear branch of the 
underlying constitutive curve is rendered unstable by a 
coupling between flow and microstructure. Within this 
model, we have found a rich variety of spatio-temporal os- 
cillatory and rheochaotic flows, many resembling experi- 
mental observations in shear banding systems. It remains 
an open challenge to delineate more fully the spectrum 
of mechanisms governing rheochaotic banded states. Ex- 
tension to higher dimensions, allowing fluctuations along 
the plane of the interface would be interesting. Within 
the present one-dimensional approach, generalisation to 
shear thickening systems [24J , curved Couette geometries 
and vorticity-banding is underway [3lj . 
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